Chaotic response of global climate to long-term solar forcing variability 
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It is shown that global climate exhibits chaotic response to solar forcing variability in a vast 
range of timescales: from annual to multi-millennium. Unlike linear systems, where periodic forcing 
leads to periodic response, nonlinear chaotic response to periodic forcing can result in exponentially 
decaying broad-band power spectrum with decay rate Te equal to the period of the forcing. It is 
shown that power spectrum of a reconstructed time series of Northern Hemisphere temperature 
anomaly for the past 2,000 years has an exponentially decaying broad-band part with Te ~ 11 yr, 
i.e. the observed decay rate Te equals the mean period of the solar activity. It is also shown that 
power spectrum of a reconstruction of atmospheric CO2 time fluctuations for the past 650,000 years, 
has an exponentially decaying broad-band part with Te — 41, 000 years, i.e. the observed decay rate 
Te equals the period of the obliquity periodic forcing. A possibility of a chaotic solar forcing of the 
climate has been also discussed. These results clarify role of solar forcing variability in long-term 
global climate dynamics (in particular in the unsolved problem of the glaciation cycles) and help in 
construction of adequate dynamic models of the global climate. 

PACS numbers: 92.60.Iv, 92.30.Np, 05.45.Tp, 05.45.Gg 
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INTRODUCTION 

Behavior of a chaotic system can be significantly 
altered by applying of a periodic forcing. Already 
pioneering studies of the effect of external periodic 
forcing on the first Lorenz model of the chaotic climate 
revealed very interesting properties of chaotic response 
(see, for instance, [H,[l|,[ll)- The forcing does not 
always have the result that one might expect 
The climate, where the chaotic behavior was discovered 
for the first time, is still one of the most challenging 
areas for the chaotic response theory. One should 
discriminate between chaotic weather (time scales up to 
several weeks) and a more long-term climate variation. 
The weather chaotic behavior usually can be directly 
related to chaotic convection, while appearance of the 
chaotic properties for the long-term climate events is 
a non-trivial and challenging phenomenon. It seems 
that such properties can play a significant role even 
for glaciation cycles, i.e. at least at multi-millennium 
time scales Cyclic forcing, due to astronomical 

modulations of the solar input, rightfully plays a central 
role in the long-term climate models. Paradoxically, it is 
a very non-trivial task to find imprints of this forcing in 
the long-term climate data. It will be shown in present 
paper that just unusual properties of chaotic response 
are the main source of this problem. 



GLOBAL TEMPERATURE RESPONSE TO 
PERIODIC SOLAR FORCING 



Figure 1 shows a reconstruction of Northern Hemi- 
sphere temperatures for the past 2,000_years (the data 
for this figure were taken from Ref. [9|). This multi- 
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FIG. 1: A reconstruction of Northern Hemisphere tempera- 
ture anomaly for the past 2,000 years. 



proxy reconstruction was performed by the authors of 
Ref. [10| using combination of low-resolution proxies 
(lake and ocean sediments) with comparatively high- 
resolution tree-ring data. Figure 2 shows a power spec- 
trum of the data set calculated using the maximum en- 
tropy method, because it provides an optimal spectral 
resolution even for small data sets. The spectrum ex- 
hibits a wide peak indicating a periodic component with 
a period around 22 y, and a broad-band part with expo- 
nential decay: 



E{f) 



(1) 



A semi-logarithmic plot was used in Fig. 2 in order to 
show the exponential decay more clearly (at this plot the 
exponential decay corresponds to a straight line). Both 
stochastic and deterministic processes can result in the 
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FIG. 2: Spectrum of the data, shown in Fig. 1, in semi- 
logarithmic scales. The straight line indicates the exponential 
decay Eq. (1). 
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FIG. 3: A sketch of the complex time plane for the parametric 
modulation with period Te. 



broad-band part of the spectrum, but the decay in the 
spectral power is different for the two cases. The ex- 
ponential decay indicates that the broad-band spectrum 
for these data arises from a deterministic rather than a 
stochastic process. For a wide class of deterministic sys- 
tems a broad-band spectrum with exponential deca y is a 
generic feature of their chaotic solutions Refs. (lH-[l4j. 

Nature of the exponential decay of the power spectra 
of the chaotic systems is still an unsolved mathematical 
problem. A progress in solution of this problem has been 
achieved by the use of the analytical continuation of the 
equations in the complex domain (see, for instance, [lit). 
In this approach the exponential decay of chaotic spec- 
trum is related to a singularity in the plane of complex 
time, which lies nearest to the real axis. Distance be- 
tween this singularity and the real axis determines the 
rate of the exponential decay. For many interesting cases 



chaotic solutions are analytic in a finite strip around the 
real time axis. This takes place, for instance for attrac- 
tors bounded in the real domain (the Lorentz attractor, 
for instance). In this case the radius of convergence of 
the Taylor series is also bounded (uniformly) at any real 
time. If parameters of the dynamical system fluctuate pe- 
riodically around their mean values with period Tg, then 
the restriction of the Taylor series convergence (at certain 
conditions) is determined by the period of the parametric 
modulation, and the width of the analytic strip around 
real time axis equals Te/27r (cf. Fig. 4). Let us consider, 
for simplicity, solution u{t) with simple poles only, and 
to define the Fourier transform as follows 

u(w) = (2^)-i/2 C''^ dt e-'^'u{t) (2) 

Then using the theorem of residues 

u(uj) = i(27r)^/^ ^ Rj e:Kp{iujXj — \ujyj\) (3) 
j 

where Rj are the poles residue and Xj + ii/j are their lo- 
cation in the relevant half plane, one obtains asymptotic 
behavior of the spectrum E{uj) — \u{uj)\'^ at large lu 

E{uj) ~ exp(-2|wy™„|) (4) 

where ymm is the imaginary part of the location of the 
pole which lies nearest to the real axis. Therefore, expo- 
nential decay rate of the broad-band part of the system 
spectrum, Eq. (1), equals the period of the parametric 
forcing. 

The chaotic spectrum provides two different character- 
istic time-scales for the system: a period corresponding to 
fundamental frequency of the system, T/„„ , and a period 
corresponding to the exponential decay rate, Tg = l//e 
(cf. Eq. (1)). The fundamental period Tfun can be esti- 
mated using position of the low- frequency peak, while the 
exponential decay rate period Tg = 1 //e can be estimated 
using the slope of the straight line of the broad-band part 
of the spectrum in the semi-logarithmic representation 
(Fig. 2). From Fig. 2 we obtain T/^ ~ 22 ± 2y and 
Tg ~ 11 ± ly (the estimated errors are statistical ones). 
Thus, the solar activity period of 11 years is really a dom- 
inating factor in the chaotic temperature fiuctuations at 
the annual time scales, although it is hidden for linear 
interpretation of the power spectrum. In the nonlinear 
interpretation the additional period Tf^n — 22y might 
correspond to the fundamental frequency of the under- 
lying nonlinear dynamical system. It is surprising that 
this period is close to the 22y period of the Sun's mag- 
netic poles polarity switching. It should be noted that 
the authors of Ref. pjj] found a persistent 22y cyclicity 
in sunspot activity, presumably related to interaction be- 
tween the 22y period of magnetic poles polarity switching 
and a relic solar (dipole) magnetic field. Therefore, one 
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FIG. 4: Spectrum of the data, shown in Fig. 1, in In-ln 
scales. The dashed straight hne indicates the Kohnogorov- 
hke spectrum: E{f) ~ f~^^^ ■ The high-frequency part has 
been cut in order to show the low- frequency part. 



sponds to so-called Kolmogorov's cascade process. This 
process is very universal for turbulent fluids and plasmas 
13], [3] ■ For turbulent processes on Earth and in Helio- 
sphere the Kolmogorov-like spectra with such large time 
scales cannot exist. Therefore, one should think about a 
Galactic origin of Kolmogorov turbulence (or turbulence- 
like processes [l^) with such large time scales. This is 
not surprising if we recall possible role of the galactic cos- 
mic rays for Earth climate (see, for instance, [2ol|-[23t). 
In order to support this point we show in figure 5 spec- 
trum of galactic cosmic ray intensity at the Earth's orbit 
(reconstruction for period 1611-2007yy [24|, cf. also Ref. 
[25|). One can compare Fig. 5 with Fig. 4 correspond- 
ing to the global temperature anomaly fluctuations. It 
should be noted, that the above discussed response of the 
global temperature to the solar activity cycles can also 
have a cosmic rays variability as a transmitting agent. 
Indeed, the interplanetary magnetic field strongly inter- 
acts with the cosmic rays. Therefore, the change in the 
interplanetary magnetic field intensity, due to the solar 
activity changes, can affect the Earth climate through 
the change of the cosmic rays intensity and composition 
(seeRefs. H-Q)- 
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FIG. 5: Spectrum of the galactic cosmic ray count rate fluc- 
tuations in the In-ln scales. The straight dashed line is drawn 
to indicate the Kolmogorov scaling law E{f) ~ j(~5/3) 



cannot rule out a possibility that the broad peak, in a 
vicinity of frequency corresponding to the 22y period, is 
a quasi-Zmear response of the global temperature to the 
weak periodic modulation by the 22y cyclicity in sunspot 
activity. I.e. strong enough periodic forcing results in 
the non-linear (chaotic) response whereas a weak peri- 
odic component of the forcing can result in a quasi-linear 
(periodic) response. 

Figure 4 shows the same power spectrum as in Fig. 2 
but in In-ln scales. A dashed straight line in this figure 
indicates a scaling with '-5/3' exponent: E{f) ~ f~^^^- 
Although, the scaling interval is short, the value of the ex- 
ponent is rather intriguing. This exponent is well known 
in the theory of fluid (plasma) turbulence and corre- 



LONG-TERM TEMPERATURE RESPONSE TO 
CHAOTIC SOLAR FORCING 

It is interesting, that the solar activity itself is chaotic 
at the multi-decadal time-scales [1^. Figure 6 shows a 
spectrum of a long-range reconstruction of the sunspot 
number fluctuations for the last 11,000 years (the data, 
used for computation of the spectrum, is available at 
[2?!). The semi- logarithmic scales and the straight line 
are used in this figure to indicate the chaotic nature of 
the sunspots number fluctuations at multi-decadal time 
scales. The slope of the straight line provides us with the 
characteristic chaotic time scale Tg ~ 176y for the chaotic 
solar activity fluctuations [2^. It should be noted that 
the 176y period is the third doubling of the fundamental 
solar period 22y, which we have discussed above (see also 
[2^ and references therein). As one can conclude from 
Fig. 4 the galactic turbulence has more profound effect 
on the global temperature than the solar forcing vari- 
ability at these time-scales, at least for the two last mil- 
lenniums. However, turbulence is a highly intermittent 
phenomenon (see, for instance, review [28|). Therefore, 
during the last two millenniums (presented in the Fig. 
4) the solar system could pass through a very intensive 
galactic turbulence patch, while before this (and after 
this) the nearest galactic environment could be rather 
quiet. This can allow us to investigate the question: 
Whether chaotic parametric modulation of a nonlinear 
system (in our case the Earth climate) also results in the 
chaotic response (in previous Chapter we have studied 
periodic parametric modulation of this system) . Figure 7 
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FIG. 6; Spectrum of the sunspots number fluctuations in the 
semi-logarithmic scales for the reconstructed data for the last 
11000 years (the data have been taken from Ref. [131 )■ The 
straight line is drawn to indicate the exponential law Eq. (1). 




FIG. 8: Spectrum of the data, shown in Fig. 7, in semi- 
logarithmic scales. The straight line indicates the exponential 
decay Eq. (1). 
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FIG. 7: A reconstruction of Antarctic temperature for the 
past 10,000 years (EPICA Dome C Ice Core ;29i]). 

shows a reconstruction of Antarctic temperature for the 
past 10,000 years using a high-resolution deuterium data 
set at EPICA Dome C Ice Core (the data for this figure 
were taken from Ref. [29]). Figure 8 shows spectrum 
corresponding to these data. The slope of the straight 
line provides us with Te — 176y (cf Fig. 6). Thus one 
can conclude that in this case we observe the chaotic re- 
sponse to the parametric modulation of the climate by 
the chaotic fluctuations of the solar activity. 



THE PROBLEM OF GLACIATION CYCLES 

The angle between Earth's rotational axis and the 
normal to the plane of its orbit (known as obliquity) 



varies periodically between 22.1 degrees and 24.5 degrees 
on about 41,000-year cycle. Such multi-millennium 
timescale changes in orientation change the amount 
of solar radiation reaching the Earth in different lat- 
itudes. In high latitudes the annual mean insolation 
(incident solar radiation) decreases with obliquity, while 
it increases in lower latitudes. Obliquity forcing effect 
is maximum at the poles and comparatively small in 
the tropics. Milankovic theory suggests that lower 
obliquity, leading to reduction in summer insolation and 
in the mean insolation in high latitudes, favors gradual 
accumulation of ice and snow leading to formation of an 
ice sheet Q. The obliquity forcing on Earth climate is 
considered as the primary driving force for the cycles of 
glaciation (see for a recent review [s^). Observations 
show that glacial changes from -1.5 to -2.5 Myr (e arly 
Pleistocene) were dominated by 41 kyr cycle (7t.[3l|.[3^. 
whereas the period from 0.8 Myr to present (late 
Pleistocene) is characterized by approximately 100 kyr 
glacial cycles [l^,!!^]- While the 41 kyr cycle of early 
Pleistocene glaciation is readily related to the 41 kyr 
period of Earth's obliquity variations the 100 kyr period 
of the glacial cycles in late Pleistocene still presents a 
serious problem. Influence of the obliquity variations on 
global climate started amplifying around 2.5 Myr, and 
became nonlinear at the late Pleistocene. Long term 
decrease in atmospheric CO2, which could result in a 
change in the internal response of the global carbon 
cycle to the obliquity forcing, has been mentioned 
as one of the principal reasons for this phenomenon 
(see, for instance, [33|-[33])- Therefore, investigation 
of the historic variability in atmospheric CO2 can be 
crucial for understanding the global climate changes at 
millennial timescales. Figure 9 shows a reconstruction 
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FIG. 9; A reconstruction of atmospheric CO2 based on deep- 
sea proxies, for the past 650kyr. 



although it is hidden for linear interpretation of the 
power spectrum. In the nonlinear interpretation the 
additional period Tfun — 100 kyr might correspond to 
the fundamental frequency of the underlying nonlinear 
dynamical system and it determines the apparent 100 
kyr 'periodicity' of th e g laciation cycles for the last 



650 kyr (cf Refs. [8| , [36| , [39| and references therein). 



And again, as in the above considered case of the 
global temperature fluctuations, one cannot rule out a 
possibility that the broad peak, in a vicinity of frequency 
corresponding to the 100 kyr period, is a quasi-Zmear 
response of the atmospheric CO2 to the weak periodic 
modulation by the 100 kyr cyclicity in the orbital 
eccentricity variations [13]. I.e. again, strong enough 
periodic forcing results in the non-linear (chaotic) 
response whereas a weak periodic component of the 
forcing can result in a quasi-linear (periodic) response. 
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FIG. 10; Spectrum of atmospheric CO2 fluctuations for the 
data shown in Fig. 9 

of atmospheric CO2 based on deep-sea proxies, for the 
past 650kyr (the data taken from [38|). Resolution of the 
data set is 2kyr. Fluctuations with time-scales less than 
2ky r could be rather large (statistically up to 308ppm 
[38j). but they are smoothed by the resolution. Figure 10 
shows a power spectrum of the data set calculated using 
the maximum entropy method. The spectrum exhibits 
a peak indicating a periodic component (the arrow in 
the Fig. 10 indicates a lOOkyr period) and a broad-band 
part with exponential decay. A semi-logarithmic plot 
was used in Fig. 10 in order to show the exponential 
decay more clearly (at this plot the exponential decay 
corresponds to a straight line). From Fig. 10 we obtain 
Tfun — 95 ± 8 kyr (the peak is quite broad due to small 
data set) and Tg ~ 41 ± 1 kyr (the estimated errors are 
statistical ones). Thus, the obliquity period of 41 kyr is 
still a dominating factor in the chaotic CO2 fluctuations. 
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